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Abstract. In this report, we discuss the measurement of the hadronic decay 
. . . modes of resonances in relativistic heavy- ion cohisions, emphasizing on RHIC 

, results. The study of resonances can provide: (1) the yield and spectra of more 

I particles with different properties to study whether the system is in thermal 

^SJ ' equilibrium; (2) an independent probe of the time evolution of the source from 

chemical to kinetic freeze-out and detailed information on hadronic interaction at 
^ ' later stage; (3) the study of medium effects at a late stage of the collisions; (4) 

. a measurement of flavor and baryon/meson dependence of particle production at 

lyQ I intermediate pt (> 2 GeV/c). 

(N ; 

(N : 

0.1. Hadron Production in Heavy- Ion Collisions 

O Historically, the discovery of several important resonances (S(1385), K*{892), 

^ . p°(770), 7], etc.) helped confirm the quark model (the Eightfold way) in the early 

2 ! 60's pp. Among more than four thousand particles discovered, most of them are 

. resonances 0. In a simple system such as e~^e~ collisions at i/i = 91 GeV, the 

(U '. production characteristics of 46 different particles have been measured With 



X 



same way of counting, only 13 particles (vr^, 7r°, K"^, K^, p°, rj, (p, K*'^, J/'^jp, A, E, Q) 
!^ ! have been measured in ultrarelativistic heavy-ion collisions up to date. 

The soft processes of QCD hadronization and interaction are not calculable with 
a perturbative approach and instead rely on phenomenological models to describe 
the particle production. In a simple system, such as e^e~ collisions, the string and 
5^ , cluster fragmentation can be used in Monte Carlo Models or in an empirical model to 
describe the data |3]. In heavy- ion collisions, sufficient interactions among particles 
after their creation may lead to thermalization and therefore a thermal model could 
be used to describe the data. In fact, thermal models have been able to successfully fit 
the data and extract thermal parameters such as temperature and chemical potential 
at AGS, SPS and RHIC energies [S]. In these models, a large fraction of the measured 
yields of stable particles comes from the feeddown of resonance decays. In some cases, 
it is the dominant contribution to the measured value. For example, in JETSET and 
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thermal models, less than 20% of measured pions (30% kaons) do not originate from 
resonances [HE], while the resonance production in elementary collisions is dominated 
by primary production. It is therefore important to measure resonances to constrain 
the feed-down and to continue testing the validity of thermal models. The same 
argument is applicable to transport models (such as UrQMD 6J, AMPT f7|) and 
hydrodynamic models. 

The first K* measurement by STAR on K''{892)/h~ ^ at ^s^= 130 GeV/c is 
consistent with thermal fit [3] with the combined statistical error and systematic error 
of the order of 30%. The ratio K*^ / K is less model dependent than K*^/h^ since 
both particles have similar quark content and differ only in their spin and mass. The 
right panel in Fig. ^ shows that this ratio is independent of beam energy in pp, pp and 
e~^e~ . Thermal model predicts this ratio to be 0.32 at chemical freeze-out at RHIC 5j. 
The STAR results on is:*(892)//s:= 0.26±0.03±0.07 in central Au+Au collisions 
at ^/s^ = 130 GeV/c and K*{892)/K = 0.20 ± 0.01 ± 0.03 in central Au-^Au 
collisions at y/s^=200 GeV/c with improved statistical and systematic uncertainty 
[8j show that this ratio does have a centrality dependence and is significantly lower 
than thermal model prediction or results from pp collisions measured by the same 
experiment. Siginificant decrease of A(1520) in Au+Au collisions with respect to pp 
and thermal calculation are also observed i.lOI. 



g 
(5 





1 1 1 1 1 1 1 1 1 1 1 1 

-B- K*° / K in AuAu 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
^<t/K*° inAuAu 


1 ' 




-■- K*° / K in pp 


/ K*" in pp 






STAR Preliminary 










1 { 




- T 

- • 








7 i 










□ 


a 




\ , 


. , 1 , , , , 1 , , , , 


... 1 .... 1 .... 1 . 





100 200 



300 400 500 600 700 

Number of Charged Hadrons 



STAR Preliminary 3) K /K ratio 



4 



A a 



Figure 1. Left panel: K*^ / K ratio and 4)/K* ratio as function of charged hadron 
multiplicity in Au+Au collisions at ^s„„ =200 GeV. The filled symbols are results 
from pp collisions. The line is a prediction from thermal model [H]. Right panel: 
K*^ / K ratio and (ji/K* ratio as function of beam energy in e+e" , pp, pp and Au+Au 
collisions. Filled symbols are for Au+Au collisions and the open symbols arc for 
elementary collisions. 



Whether the remarkable success of the statistical model may simply mean phase- 
space dominance and the "temperature" and "chemical" potentials are nothing but 
Lagrange multiplier characterizing the phase-space integral JT] is still unclear. The 
existence of flow is considered a strong evidence that certain amount of rescattering 
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takes place among the particles in the bulk matter formed in heavy-ion collisions. 
However, as stated in Ref . [TTj , "to which extent they are sufficient to form matter in 
the Boltzmann sense is, however, not clear" . The strong rescattering and regeneration 
of resonances after chemical freeze-out as discussed are direct evidence of rescattering 
of particles in the bulk matter formed in heavy-ion collisions. It is difficult to imagine 
that such interaction will not be stronger at even earlier stage. 

0.2. Time Evolution of the System 

Resonances which decay into strongly interacting hadrons inside the dense matter 
are less likely to be reconstructed due to the rescattering of the daughter particles. 
On the other hand, resonances with higher pt have a larger probability of decaying 
outside the system and, therefore, are more likely to be reconstructed. Examples of 
how these measurements can be used as a signature of freeze-out djTiamics are shown 
in Refs. [Ellin!. 

Recent transport model (UrQMD) calculation shows significant modification 
of resonance population and pr spectrum at mid-rapidity due to rescatterings [H]. 
In addition, the resonance yield could be increased during the rescattering phase 
between chemical freeze-out (vanishing inelastic collisions for both stable particles and 
resonances) and kinetic freeze-out (vanishing elastic collisions) via, for example, the 
elastic process nK K*^ — > ttK. This regeneration mechanism partially compensates 
for resonance decays during a possible long expansion of the system and increases the 
observable ratio K*^ / K. The tttt scattering is the dominant process in destroying the 
reconstruction of at later stage while the Kn reaction is the dominant channel 
of i^* regeneration. The tttt cross section is larger than the Kt^ cross section by a 
factor of ~5 jHj around the resonances. This means that at sufficiently late stage, 
the rescattering of the daughters via the scattering of decay pions will be much larger 
than the regeneration process via the Ktt reaction. Therefore, significant effect from 
these rescattering should be detectable experimentally. 

It has been shown in this conference that there is a significant decrease of 
K* /K pj and A*/A pH] ratios while (p/K ratio is independent of beam specie and 
centrality This indicates that the yields are related to the cross section of 

the interactions between their decay daughters and the surrounding hadrons. It can 
also be seen from the px spectra of and as shown in right panel of Fig. IHl that 
the low Pt part of is suppressed relative to and other particles. Resonances 
being discussed are 0(1020), fs:*±(892), A(1520), E(1385), A++(1232), p°(770) and 
/o(980). The ordering of the lifetime is 0(1020) > A(1520) > S(1385) > K*{SQ2) > 
/o(980) > A++(1232) > p°(770). The ordering of their regeneration cross section 
is A++(1232) > p°(770) > /o(980) > K*{892) > S(1385) > A(1520) > 0(1020). 
Meanwhile, the rescattering of daughters in all cases is dominated by either tttt or pn. 
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The ordering of their mass difference to ground state is S(1385)/A< A++(1232)/p< 
K*{892)/K < A(1520)/A < p^/vr < /q/tt. If the rescattering is the dominant 
process, the resonance yields and pr spectra will be related to the lifetime and cross 
sections. On the other hand, if the dominant processes are both the regeneration 
and rescattering with a long time expansion which tend to equilibrate the yields at 
kinetic freeze-out, the resonance yields and pt spectra will be related to the Boltzmann 
factor. By systematically comparing the yields and px distributions of resonances with 
other particles it may be possible to distinguish different freeze-out conditions; e.g. 
sudden freeze-out [T^ IT^ or a slow hadronic expansion |17j . From the preliminary 
results of / R- , A{1520) /A, /n^ Jo/n- and (p/K' [IHl El CHI EE] measured by 
STAR and the comparison with thermal model predictions at chemical freeze-out 
jHj, it seems that the ratios follow the regeneration cross section. This implies 
that large rescattering of daughters happens after chemical freeze-out. However, the 
regeneration is not enough to compensate the loss to equilibrate the yield due to fast 
freeze-out and/or decreasing temperature. 

It has been noted that jHJ [1211201 due to the overpopulation of tt's at kinetic freeze- 
out, the resonances having a tt daughter will have higher yield due to additional pion 
chemical potential. This will change the ordering of the particle ratios mentioned 
above. With these precise measurements, it is now clear that the yield of resonances 
are not consistent with thermal predictions E] at chemical freeze-out. The 
interpretations of these results are being actively pursued as seen in these proceedings 
and recent publications. The results from the combination of a thermal model as 
the initial condition together with the UrQMD transport model compared to the 
measured yields and spectra of resonances and other particles will constrain both the 
thermal model and evolution of the system [HI . Additional results from A"'"^ and 
S(1385) are coming soon from STAR as well STAR is also carrying out the 

studies of resonance elliptic flow {v2) and decay angular distribution (cos 6'*) which 
are sensitive to the geometry of the system and the rescattering [THIini- These studies, 
however, require much more statistics and better systematical error than the yields 
and pt spectra. 

0.3. Medium Effect in Hot and Dense Matter 

Among the proposed signals of phase transitions in hot, dense nuclear matter 
produced via relativistic heavy-ion collisions are modifications of the meson resonance 
production rates and their in-medium properties fHl- When the resonance lifetime 
is comparable to the evolution time scales of the phase transition, the measured 
properties associated with the resonance (such as mass, width, branching ratio, yield 
and transverse momentum (pt) spectra) will depend upon collision dynamics and 
chiral properties of the medium at the high temperature and high energy density. 
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The famous result is the enhancement of dileptons in the medium mass range around 
500 MeV by CERES/SPS Models with in- medium mass modification of 

were proposed to explain the large enhancement. The resonances measured in the 
hadronic decay channel probe only the late stages of the collision. It has been argued 
that significant effect on the reconstructed resonance mass can exist in such late stage 
due to phase space [201 1211 , interference [23] , rescattering [Bj and dynamical effects 
|20| I19j. As shown in ^§110] and here in Fig. |21 STAR has observed significant mass 
shifts in both p° — > vt+tt" and K* — > ttK decays [HI ^]. The mass positions are 
consistently lower than the PDG (Particle Data Group [H]) values and pt dependent 
with a large discrepancy at lower pt and with the measured values approaching the 
PDG values at higher pT. It was found that the p° mass at low p^ was about -30 
MeV/c^ and -70 MeV/c^ below the PDG values in pp and peripheral Au-|-Au collisions 
while the K*^ value is about -15 MeV/c^ below the PDG values. The statistical error 
of K*^ in Au-|-Au collisions is too large to be conclusive. The mass and width of other 
resonances (0, E* and A*) and strange particles {Kg, A) are consistent with the PDG 
values within the STAR experimental resolution [THl I10| 1^. 
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Figure 2. Left panel: mass as function of pt in Au+Au collisions with centrality 
of 40-80%. The filled circles are results from the fit to a Breit-Wigner function 
modified by a phase-space factor. Right panel: if*" mass vs pT in pp and Au+Au 
collisions. 



The first question one should ask is whether the mass shift of these particles 
(especially the copiously produced p° vector meson) has been observed anywhere. In 
fact, previous measurements of the p° meson in hadronic Z° decays by the OPAL, 
DELPHI and ALEPH collaborations indicate that the p° line shape is considerably 
distorted from a relativistic p-wave Breit-Wigner shape, especially at relatively low 
momentum in multipion systems [23 l2Sll2niin2I • A mass shift that may be -30 MeV/c^ 
or larger was observed [2111211 • The OPAL experiment at LEP has reported mass 
shifts of -10 to -30 MeV/c^ in the position of the maximum of the resonance consistent 
with the observed p° mass shift [23 • In the paper by the DELPHI collaboration, the 
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low momentum region, below 2.3 GeV/c was excluded and a Breit-Wigner shape 
was used to fit the p° |28| Although the fitted mass was some five standard 
deviations below the PDG value, rates for the p° were nonetheless extracted using a 
conventional Breit-Wigner shape fIE\ • The DELPHI collaboration reported the 
mass as 757±2 MeV/c^ j2H], which should be compared to the PDG average value of 
775.9 ±0.5 MeV/c^ from r decays and e~^e~ interactions. The ALEPH collaboration 
also reported that the p° resonance in the data seemed to be shifted towards lower 
masses pUj . 

In pp collisions, the meson has been measured at ^/s = 27.5 GeV jSU]. This 
is the only pp measurement used in the hadroproduced p° mass average reported in 
the PDG. In this measurement, a mass of 0.7626 ± 0.0026 GeV/c^ was obtained 
from a fit by a relativistic p-wave Breit-Wigner function times the phase space piOJ| . 
However, if the p^ is fit only to a relativistic p-wave Breit-Wigner function, a mass 
shift of -30 MeV/c^ is observed. 

In what follows, I would like to discuss some of the predicted effects on the mass 
and width of the p" in the medium. 

0.3.1. Phase Space and Rescattering If the resonance observed is generated from 
scattering (e.g. tt+tt^ p. Kit K*), then the mass distribution will have 
contribution from the phase space population of the parents (which are also 
daughters). This distribution is characterized by the Boltzmann distribution in a 
thermal system that tends to populate the low invariant mass of the resonances and 
therefore effectively shifts the resonance peak to lower values. Since in the invariant 
mass distribution (without contamination from particle misidentification often present 
in the experiments) the parents (daughters) distribution represents the phase space 
and the experimental acceptance and efficiency, the p° mass is often obtained by 
fitting the same event distribution of tt+tt^ to 

BG + PS X BW = BG + BGx BW = BG{1 + BW), (1) 

where BG corresponds to the background, PS is the phase space factor and BW is the 
relativistic p-wave Breit-Wigner (e.g. Ref. [2111231201 which was referred to in PDG). 
However, this has a strong assumption that all resonances are created by final state 
rescattering. In elementary collisions {e~^e~,pp and pp), some of the resonances are 
directly from string fragmentation which should not be modified by the phase space 
population of other particles. In heavy-ion collisions, such a fit is better justified since 
the system is much larger and lives longer. This thermal factor has been discussed 
in several papers [2111 El]- Usually, this constribution to the total mass shift is quite 
small (<20 MeV) and scales roughly with P^. 

A more detailed study can be done using transport models, such as UrQMD 6J. 
The left panel in Fig. |3] shows the p° mass for different centrality from UrQMD in 
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Au+Au collisions. The mass shift is about -30 MeV/c^ for central Au+Au collisions 
and vanishes for peripheral collisions. 

0.3.2. Interference and Bose-Einstein Effect There are several sources of t!:^ti~ 
correlations around the p° mass: 

• Direct p° production and its decay via ^7r"'"7r~; 

• TTTT scattering via p° resonant state; tt+tt^ ^p^7r+7r~; 

• Direct tttt scattering via tt+tt^ ^tt+tt^; 

• Interference between previous two scatterings; 

• Other TTTT waves (S-, P-, D- waves). 

The interference has been studied in great detail in the coherent p° photoproduction. 
The STAR Collaboration used the same detectors and the same beam to trigger on 
the p° photoproduced in ultraperipheral Au+Au collisions. The mass and width 
are consistent with the PDG values and the interference strength is consistent with 
other experiments jSIj. It has been argued that a similar interference may happen in 
p° production in tt^tt~ scattering This, however, may not be compatible with 
vector moson dominance (VMD) |34j. The p" line shape in multihadronic decay is 
addressed in details in Ref. The final state Bose-Einstein effect was introduced 
to explain the p° hue shape measured at LEP [23 EHl EHl IS3 CHI • 

0.3.3. p" Mass In-Medium Modification Although the effects mentioned above are 
interesting effects on their own, in heavy-ion collisions the most interesting effect 
is to identify the modification of particle (vector meson) properties in the medium. 
The modification of the hadron spectral distribution even at dilute systems are of 
particular interest since it will be a precursor toward the chiral phase transition [.19j. 
A couple of calculations using Brown-Rho scaling or dynamic width broadening can 
qualitatively explain the mass shift fOl 1201 • The preliminary results from STAR 
clearly should be compared with theoretical calculations more quantitatively in order 
to understand to what extent the measurements reflect a true in-medium modification 
[T^ . The hadronic channel measured by STAR coupled with future measurements 
of the dilepton channel of vector mesons will provide an unique tool to study the 
in-medium effects as well as the properties of the hot and dense matter. 

0.4. Mass or Baryon/Meson Dependence at Intermediate pt 

A very prominent effect seen at RHIC is that the pt spectra of different particles 
are very different from that of pp collisions at intermediate pt f33.,. The effect has 
been attributed to either a baryon-meson effect or a mass dependence from radial 
flow. Unfortunately, the baryons are often heavier than mesons (protons are heavier 
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than pions and A's are heavier than kaons). The meson resonances such as (p and 
K* are as heavy as protons, but have the properties of mesons (two valence quarks). 
Therefore, it is important to study the pt spectra of these particles to distinguish 
baryon-meson or mass effect. The left panel in Fig. El depicts the particle spectra of 
7r,K,p,A,K*°,(f) and S measured by all four RHIC experiments [SHI CHI El CEl ESI EEl ■ 
The right panel in Fig. El shows the pt spectra of particles with similar mass in 
central Au+Au collisions measured by STAR. The p spectrum includes the hyperon 
decay feed-down. The spectra of p, and S are scaled by a factor of 0.5, 2.5 and 5 
respectively. It is surprising to see that these spectra are very close to each other in 
the intermediate pt range between 1 GeV/c and 4 GeV/c. Although the rescattering 
at late stages suppresses the K* yield, the intermediate pt region {prD--^ GeV/c) 
has a smaller effect [Hj. For pt< 1-0 GeV/c, we can see that particles with different 
properties deviate from each other. For example, has a smaller hadron cross section 
and therefore has steeper slope at low px, while the i^*" daughters rescatter in the 
hadronic matter and the low pt K*^^s are suppressed. A study of the (p/K* ratio 
at intermediate pt will reveal the production mechanism, e.g. fragmentation and 
recombination at partonic level |HZ|- From Fig. El we see that the (p/K*^ is about 
0.4. This is similar to the value in pp collisions, as shown in Fig. ^ This indicates 
that the changes from pp to central Au+Au collisions are similar for the and K*^. 
Further comparison between experimental data and theoretical predictions on these 
spectra are needed. 




(GeV/c) 



Figure 3. tt, K, p, A, K*'^, (f> and ^ pr spectra in central Au+Au collisions measured 
at RHIC. See text for details. 

Since there are more quark than antiquark jets and final hadrons from gluon 
jets should have equal number of particles and antiparticles, one should find more 
p. A, K and K* than p. A, K and K* at mid-rapidity in heavy-ion collisions. It was 
found that the flavor dependence and suppression factor would be a good probe of 
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the energy loss and other novel effects tHB] • Fig. |3] shows the HIJING simulation 
of antiparticle to particle ratio as function of pt at midrapidity {\y\ <0.5) in Au+Au 
collisions at y^s^ = 200 GeV. It can be seen that the ratios of K/K and K*/K* 
have the same sensitivity to the flavor dependence of particle production as p/p 
and A/A. In addition, the K*^ measurement will allow the measurement of the 
isospin effect by using K*^/K* ~ 1.0 and comparing to the flavor dependence of 
K*/K*. Experimentally, the momentum resolution of lower px particles is better 
than higher px due to the rigidity (curvature) in the magnetic field (e.g. STAR). The 
reconstruction of two or three lower momentum daughters to form a high p^ particle 
is less sensitive to the distortion when we want to study the flavor and isospin effect. 
Experiments have shown that the ratios of K^/K^ and A/A at low px do not 
depend on px- It is important to find the onset (on px) of the strong fiavor dependence 
as predicted. 




Figure 4. Left: p'^ mass integrated over all px from UrQMD vs centrality in Au+Au 
collisions. Right: Effect of valence quark in particle ratio. HIJING simulation of 
antiparticle to particle ratio as function of pr- In case of K*, we have K*/K* 
(flavor) and K*+/K*° (isospin). 



In summary, studies of strongly interacting resonant states open new approaches 
to the study of relativistic heavy-ion collisions. We found that the yields of resonances 
deviate from thermal fits at chemical freeze-out. This is likely related to the evolution 
of the system and continuous interactions after chemical freeze-out. The observation 
of p° and K*'^ mass shifts in both pp and Au-|-Au collisions are actively investigated 
theoretically and experimentally. These results show clear dynamic effects at the 
freeze-out. The transverse momentum spectra of resonances at intermediate and 
high px can help disentangle different effects. 

The auther would like to thank the useful discussions with Drs. G. Brown, M. 
Bleicher, P. Fachini, L. Gaudichet, T. Hallman, H. Huang, R. Longacre, J. Ma, R. 
Majka, C. Markert, R. Rapp, F. Retiere, J. Sandweiss, E. Shuryak, A. Tai, T. Ullrich, 
N. Xu, H. Zhang. The references from LEP and NA27 were done together with P. 
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BNL) and P. Fachini. 

[I] L. W. Alvarez, Nobel Prize Lecture, 1968. 

[2] J. Rafelski, Quark Matter 2002, Nantes, France. 
[3] D.E. Groom et ai, Eur. Phys. J. C 15, 1 (2000). 

[4] Y. J. Pei, Z. Phys. C 72, 39 (1996); P.V. Chliapnikov, Phys. Lett. B 470, 263 (1999). 

[5] N. Xu and M. Kaneta, Nucl. Phys. A 698, 306c (2002); F. Becattini, Z. Phys. C 69, 485 (1996); 

F. Becattini and U. Heinz, Z. Phys. C 76 269 (1997); P. Braun-Munzinger et ai, Phys. Lett. 

B 518 41 (2001); J. Stachel, these proceedings; D. Magestro, CIPANP 2003. 
[6] M. Bleicher, J. Aichehn, Phys. Lett. B 530, 81 (2002); M. Behcher, these proceedings. 
[7] B. Zhang et ai, Phys. Rev. C 61 067901 (2000). 

[8] C. Adler et al, STAR CoUaboration, Phys. Rev. C 66, 061901(R) (2002). 
[9] H. Zhang for the STAR Collaboration, these proceedings. 
[10] L. Gaudichet for the STAR Collaboration, these proceedings. 

[II] V. Koch, Nucl. Phys. A 715, 108c (2003). 

[12] J. Letessier, et ai, nucl-th/0011048; G. Torrieri, J. Rafelski, Phys. Lett. B 509 (2001) 239. 
[13] W. Broniowski, W. Florkowski, B. Hiller, nucl-th/0306034, W. Broniowski, W. Florkowski; 

Phys. Rev. Lett. 87, 272302 (2001). 
[14] S.A. Bass et ai, Prog. Part. Nucl. Phys. 41 225 (1998). 

[15] C. Adler et ai, STAR Collaboration, Phys. Rev. C 65, 041901(R) (2002); E. Yamamoto for the 

STAR Collaboration, Nucl. Phys. A 715, 466c (2003); 
[16] J. Ma for the STAR Collaboration, these proceedings. 
[17] L.V. Bravina, et aL, hep-ph/0010172, D. Teaney, et ai, nucl-th/0110037| 
[18] P. Fachini, STAR Collaboration, Nucl. Phys. A 715, 462c (2003). 

[19] R. Rapp and J. Wambach, Adv. Nucl. Phys. 25 1 (2001); R. Rapp, Phys. Rev. C 63 054907 

(2001); R. Rapp, e-Print Archive:hep-ph/0305011 
[20] E.V. Shuryak, G.E. Brown, Nucl. Phys. A 717, 322 (2003). 
[21] G. Torrieri and J. Rafelski, these proceedings. 
[22] An Tai, Sevil Salur, private communications. 
[23] CERES/SPS Collaboration, Phys. Rev. Lett. 75, 1272 (1995). 
[24] P.P. Kolb and M. Prakash, Phys. Rev. C 67, 044902 (2003). 
[25] R.S. Longacre, e-Print Archive: |nucl-th /0303068 
[26] H. Long for the STAR Collaboration, these proceedings. 
[27] The OPAL Collaboration, Z. Phys. C 56, 521 (1992); EPJ 5, 411 (1998). 
[28] The DELPHI Collaboration, Phys. Lett. B 298, 236 (1993). 
[29] The ALEPH Collaboration, Z. Phys. C 69, 379 (1996). 
[30] The NA27 Collaboration, Z. Phys. C 50, 405 (1991). 
[31] The STAR Collaboration, Phys. Rev. Lett. 89, 272302 (2002). 

[32] G.D. Lafferty, Z. Phys. C 60, 659 (1996). 

[33] L Vitev, M. Gyulassy, Quark Matter 2002, e-Print Archive: |hep-ph/0208108 

[34] M. Harada and C. Sasaki, hep-ph/0306277| 

[35] T. Ullrich, Nucl. Phys. A 715, 399 (2003). 

[36] J. Castilo for the STAR Collaboration, these proceedings. 

[37] R.J. Fries et aL, nucl-th/0306027, 

[38] X.N. Wang, Phys. Rev. C 58, 2321 (1998). 

[39] C. Adler et al., STAR Collaboration, Phys. Rev. Lett. 89, 092301 (2002). 



